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ELECTROMAGNETIC PROPERTIES OF A CIRCULAR APERTURE IN A 
DIELECTRIC-COVERED OR UNCOVERED GROUND PLANE 


By M. C. Bailey, S. N. Samaddar, * and Calvin T. Swift 
Langley Research Center 

SUMMARY 

The problem of an aperture fed by a circular waveguide excited in the dominant 
transverse electric mode (TEjjj is considered. The input admittance in the aperture 
reference plane for the case of a homogeneous lossy dielectric cover of finite thickness 
and the radiation patterns for the case of no cover are computed from integral transforms, 
for which a TE^i mode variation of the electric field is assumed at the aperture. Com- 
putations for the uncoated antenna show that the patterns and aperture admittance are 
similar to those of a square aperture with a field distribution given by the TEqj rectan- 
gular mode. 

It is noted that the uncoated aperture is well matched above midband of the circular 
waveguide in contrast to the open-end rectangular waveguide. This inherent property 
would be advantageous for possible use of the open-end circular waveguide as an antenna. 

INTRODUCTION 

The antenna under consideration consists of a circular waveguide, opening onto a 
flat dielectric -covered ground plane as shown in figure 1. The dominant mode is incident 
upon the aperture, and reflected modes of higher order are assumed to be negligible. The 
fields inside the waveguide are therefore uniquely defined by a discrete TE (transverse 
electric) vector potential containing one unknown, the reflection coefficient r. 

The fields inside the dielectric slab are represented by standing waves and are 
described by integral transforms of TE and TM (transverse magnetic) potentials, each 
containing two unknown coefficients and each potential exterior to the dielectric slab con- 
taining one unknown coefficient. The boundary conditions at the aperture and at the outer 
surface of the dielectric slab are used to solve for the unknowns. This solution results 
in an integral expression for T, or equivalently, the input admittance at the aperture 

*Raytheon Company, Space and Information Systems Division. (Dr. Samaddar 's 
collaboration with the Langley authors resulted from recognition during their separate 
endeavors of mutual interest in the same problem.) 



reference plane, z = 0. (Further mention of input admittance is referred to this plane.) 
When the permittivity of the dielectric medium ej equals the permittivity of free space 
6 q, or when the thickness of the dielectric coating d is zero, the expression correctly 
reduces to that derived by Mishustin (ref. 1), who used a Green's function formalism; 
however, the integral transform method appears to be more easily applied to the slab 
problem. 

The radiation fields for the case of no dielectric coating are derived by using the 
method of stationary phase, and they are in agreement with the patterns computed by 
modifying the expressions given by Risser (ref. 2). Numerical values of patterns and 
admittance for the aperture with no coating are plotted and compared with those of a 
square aperture of equal area, excited in the TEqj rectangular waveguide mode. 

SYMBOLS 


A 2 z -component of vector potential, electric type 

z -component of vector potential, magnetic type 
a waveguide radius 

bj n normalized input aperture susceptance 

d thickness of dielectric coating 

E electric field intensity 

g in normalized input aperture conductance 

g s normalized surface wave conductance 

H magnetic field intensity 

J -F 


JiUp) 

k 11 

K z 


Bessel function of first kind of order 1 


axial TEjj wave number in the guide, 
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axial wave number in the dielectric medium 


axial wave number exterior to the structure 


propagation constant in free space, 2 tt /\ q 

1^1 

propagation constant in the dielectric medium, ^-l/— 

A o y e o 

index of refraction of dielectric slab, kj /k Q 


complex power 


r, e,<p 


spherical coordinates 


unit vector in z direction 


Cartesian coordinates 


first root of Jj' ^x^) = 0 


normalized input admittance of aperture, g in - jb ln 


admittance of free space 


normalized radial wave number, £/k 0 


location of surface wave poles in complex 0-plane 


reflection coefficient 


differential vector operator 


6(£-£) Dirac delta function 


permittivity of free space 
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e l 

? 


4o 

n 

p,<p,z 

0 ) 


complex permittivity of the dielectric medium 

radial wave number exterior to waveguide (used as a continuous dummy 
variable) 

wavelength in free space 

permeability of free space 

radial wave number exterior to waveguide 

cylindrical coordinates 

angular frequency 


Intermediate functions: 


coefficients of magnetic vector potential in the dielectric slab 
coefficients of electric vector potential in the dielectric slab 
Ti(£) coefficient of magnetic vector potential external to the structure 


t 2 («) 


coefficient of electric vector potential external to the structure 


Xj.{£;),X2(£)> x 3(£)A expressions given by equations (14) to (18) to simplify the 

X4.(l;),W(£) J solution to the boundary -value problem 

Subscripts: 


vector components in cylindrical coordinates 
r vector components in spherical coordinates 
n order of surface wave mode 

Superscripts: 

I denotes field components and vector potentials inside waveguide 
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u 


denotes field components and vector potentials in the dielectric region 


III denotes field components and vector potentials in the region outside the 

dielectric slab 

* complex conjugate, except when used in connection with the magnetic vector 

potentials 


THEORY 


Statement of Problem and Method of Approach 


The geometry of the problem is shown in figure 1. A circular waveguide, fed in 
the dominant TEji mode opens onto an infinite, perfectly conducting flat ground plane 
coated with dielectric material of thickness d and complex index of refraction N. 


The medium inside the waveguide and 
external to the structure is assumed to 
be free space. It is also assumed that 
no higher -order waveguide modes are 
excited at the aperture and that the 
fields everywhere vary in time as e"^*. 

In order to derive the input admit- 
tance expression for the circular wave- 
guide aperture, the fields in each region 
are constructed from vector potentials 
which are solutions to the wave equation 
in cylindrical coordinates. Implicit 
relations between the unknown coeffi- 
cients are obtained by equating tangen- 
tial E and H fields at the boundary, 
z = d, and requiring that tangential E 
be continuous at the aperture, z = 0. 

The integrals are evaluated with the use 
of the Dirac delta function 6 for Hankel 
transforms. The input admittance is 
then found by equating the complex 
power expressions at the aperture, and 
the solution is expressed in a form more 
convenient for numerical integration. 


y 



Figure 1 .- Circular aperture coated with 
homogeneous lossy dielectric. 
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Fields and Vector Potentials 


■Hill 


Since only the dominant TEj j mode is assumed to exist inside the waveguide, the 
fields for z = 0 can be uniquely described in terms of a single vector potential as 
follows: 



j 


whereas, for z > 0, the fields must be described by both TE and TM potentials such that 


and 



where Ag and A^, each of which satisfies the wave equation, are defined by Stratton 
(ref. 3). 

Since all higher -order TE and TM reflected modes are assumed negligible, the 
solution for z < 0 is given by 




jk 2 U z 


+ Te 


-jk z U z 


') 




fx n p 


cos 


(4) 


6 


where the argument of the radial dependence of equation (4) follows from the requirement 
that tangential E = 0 at p = a. The general solution to the wave equation in regions n 
and m requires that both sets of potentials contain a Fourier series in <p. However, 
from the boundary conditions, all terms except the coefficient of cos d> in the series for 

A *n A *m 

— — and --f — and the coefficient of sin d> in the series for and A* 111 must 

vanish. 

With this a priori information, a straightforward separation of variables and an 
integration over radial modes § leads to solutions of the form 

^-= cos $ J* + R 1 (|)e" ;|kziZ ]d4 


A z n = sin <p f° IJiUp)^?)® 315212 + R2U)e" jkz17 ] 


(5) 




and 


A *m 

Az 


= cos $ f ITj^JjUp)^ 20 ^ 

*^0 


= sin <f> |T 2 (?)J 1 (?p)e 3kz ° Z d4 


(6) 


where the Bessel function of the first kind is required to assure finite fields at p = 0. 
Equations (5) and (6) can also be obtained by performing the Hankel transform of the wave 
equations satisfied by A z and A z /e with respect to the kernel Jj(£p). 

The radiation condition is guaranteed if 


^zo 

kzo 

where k 0 2 = w 2 p 0 e 0 and 

^zl 

^zl 

where kj 2 « a> 2 M 0 €i- 


fko 2 - I 2 

(ko 2 >« 2 )] 

i(/| 2 - *o 2 

(ko 2 < £ 2 ) J 


(ki 2 > { 2 )1 

j^-kj 2 

(k! 2 < 5 2 ) J 
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When equations (4), (5), and (6) are substituted into equations (1), (2), and (3), 
respectively, explicit expressions for the fields in terms of the unknown coefficients are 
obtained. These field equations are given in the appendix. 

Evaluation of Coefficients 

The unknown coefficients now remain to be solved. For this purpose, appropriate 
boundary conditions are applied at z = 0 and z = d. Continuity of the tangential com- 
ponents of the fields at z = d leads to the four equations: 

(1) Continuity of E^: 

e fcld [ll<«)$ ? Jl(SP) - 12(0^- hM 

+ J,<S0 ♦ K 2 <«)^ ^ ^ - T 2<{>^ f J l lSp >] (9a) 

(2) Continuity of Ep: 

e K^^£l. hU) 2 **-$„(*, 

+ + It " p Jkz ° d - t 2 ( S )j^ a (9b) 

(3) Continuity of 

- a‘ ikzld [*id>!£ ^ + £ J <H ■ elkzod [ T id)^ ^ - t 2<*>£ & JlC«P)J (20 

{4) Continuity of Epi 

- b- ikzl fid>!£ £ W ♦ «2d)i hd£>] . ± im - t 2(5 ,a iiM] (9d , 
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Solving equations (9a) through (9d) by the method of determinants yields 


RlU) - + Rl<£)j < 10 > 

- ^(«)] (11 > 

Tl ( S) ' + R l<«>] < 12 > 

T 2 (4) - < 13 > 

where 

W(|) = ej(e jkzld *■ e' jkzld ) » to(e ,kzld - e' jkzld )" 

- u 2 ♦.n*,%** ld -■** ld ) < ld > 


x l(?) ~ ( k zl - k 20 ) e ^ zl ( e l k zo “ + ^l^zo * % k zJ ( 15 ) 

Xg(^) = (*1 “ e i k zo) e ^ k zl _ ^zo) + k zl * k zoj (16) 

X 3 (4) = ♦ *^ d f ♦ - e' ikkl f 

•(<A.«.4 !feld - e ' !fell l < 17 > 

X4U) = fo k zl ei kz l d e- lk ™ d k zl (e> kz l d +C -J k ^ d ) 2 - c-fel d f 

- (*zo + kzl)(e 21kzld - e‘ 2ikzld )i (18) 
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Continuity of E ^ at z = 0 gives 


A(1 + r % J i(T) * So d{ { { Il<{) + Rl<5 ^ Jl(iP) ' ~ ' X!p Zl Ji ~ 


and continuity of Ep at z = 0 gives 

.(h£) 


A( 1 + r)J A^ - ss d( ^ 


(19) 


(20) 


The coefficients I^(£) and Rj(|) may be eliminated by first multiplying equa- 
tion (19) by Ji(?p) dp and equation (20) by p dp, where £ is introduced as a 

dummy variable. Then, integrating each equation with respect to p and adding them 
gives 


r- 


A(1 + r) Jo k 




dp 


-r 

Jo 


d| 


|kzifl 2 (^) - RaU)] 


£i> 


Mo € l 


r 

Jo 


r. 


j^ Wittp) djjdp) dj x (?p) 


dp 


dp 


P dp 


( 21 ) 


where the integral over p on the right-hand side of equation (21) was extended from 
symbol a to ® because of the boundary condition, which requires tangential E = 0 
on the ground plane. 

The Integration of the left-hand side of equation (21) is trivial. However, a few 
manipulations are necessary to integrate the right-hand side. First, Bessel's equation 
is used in the following form: 


P 


d [ dJi (|P)~ 
dpj: dp 


+ (t- 2 p 2 - l)ji<?p) = 0 


(22) 


J,(£p) 

Equation (22) is then multiplied by ~~ — dp and integrated from 0 to <*>. Second 
derivatives are eliminated by integrating by parts, and the following relationship is 
obtained: 
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r 

J o 


J l(Sp)Jl(gp) + dJi(lp) djj(^p) 


dp 


dp 


p dp = i 2 r°° pjjdpjjj^pjdp 

J 0 

* £5(| - 5) (23) 


where the definition of the Dirac delta function for HanJsel transforms is given elsewhere 
(ref. 4). 

With the use of equation (23), it is easy to show that the coefficients ^(l), Rgd), 
and F are implicitly given by 

[l2tl> - Ra<«3 = " A(1 + r)cri/i o e l j 1 (x 11 )j 1 ( ?a ) (24) 

I 'kzl 

djf(fp) 

When equation (19) is multiplied by p — dp and equation (20) by Jj(JJp)dp 

and both equations are added and integrated with respect to p, an implicit expression 
for I ]_(£), Rj(|), and T is obtained in a manner analogous to the previous development 

of equation (24). The result is 


[ii(u + mail - A( n * r) i xil ^i 

- .* 




(25) 


where the prime on the Bessel function denotes differentiation with respect to the 
argument. 

Continuity of complex power across the aperture completes the boundary conditions. 
The complex power expression in region H, at the aperture, is found by integration of the 
Poynting vector over the aperture: 

P * a - i fj S~ OP *0 P ( E p n H* - *;“*?) <*» 

Substitution of the field equations from the appendix, integrating over <p, and combining 
terms yields 
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■ + + w] 

+ ^*H 2 * a) ‘ R » ( # l(c) + - wr|C^ Jl<{p)Jl(?p ^) 


P dp 


The first integral over p gives | times the delta function, 5(£ - 0 > and the 
second integral over p vanishes; hence, when equations (24) and (25) are used, equa- 
tion (27) reduces to 






[Jiteajpjwte) + 2X 2 <$j| 
W(|) “ 


(28) 


The complex power in region I at z = 0 is 


P«-1C 

2 J 0 


2jf oa 


j** dp dc^E* 1 !^ - 


F *ItJ 


a 



(i + r*)(i - r) 



(UP 0 


0 



(29) 
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which integrates to 

-n*I _ JTUI 2 


f|A| 2 (i + r*)(i - r)^jf -(gi) 2 |c n ) 2 - 


(30) 


Input Admittance 

The input admittance, normalized to the waveguide characteristic admittance, is 
found by setting P** » P* 11 at z = 0. This substitution immediately gives 

^in 3 Sin " j^in 

_ i - r 
l + f 


00 r 


! ( x n) ( 


a 2 ko 


J 




2oj 2 M 0 ei 



- 2X t U>] 


W(5) 


/ 

>2 


■ -^ 2 

\ 

' 




J^iUa)] 2 ^) + 2X 2 (0]1 

«kzlW(5) 


(31) 


Equation (31) can be simplified by examining the expressions for Xjd), XgU), 
and W(S) (eqs, (14) to (16)). This simplification, after being factored, gives; 

W(0) - 2XiO) _ ^l _ - _ 0 2 cos(k o d^N 2 - /3 2 ) - j /n 2 - /3 2 sin(k Q d^N 2 - 0 2 J 
® ^27^2 cos(k 0 d ^N 2 - j3 2 ) - j ([l - 0 2 s in(k 0 d ^N 2 -£ 2 ) 


(32) 


and 


W(ff) + 2X 2 Q3) |/n 2 7^2 cos(k 0 d^N 2 -7 2 ) - jN 2 ^l - /3 2 sin(k 0 d|/N 2 - 0 2 ) 

W ^ ) N 2 yT - ^ 2 cos(k 0 d^N 2 - /3 2 ) - j^N 2 - j3 2 sin(k 0 d^N 2 - ^ 2 ) 

where k 0 and N = kj/k 0 . 


( 33 ) 
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The input admittance expression then becomes 



(34) 

Equation (34) is basically a contour integral along the real axis of the complex 
/3-plane. For the case of a lossless dielectric coating, the singularities of the integrands 
in this equation occur on the real axis; therefore, the path of integration for numerical 
computation must be deformed around these singularities by semicircular excursions in 
the complex 0- plane. These poles are of first order and occur at the values of /3 = 0 n , 
which are the roots of the transcendental equations 

For TE modes: 

j fl - 0~ 2 = |p - 0 n 2 cot(k 0 d '/n 2 - 0^) {3 5a) 


For TM modes: 



(35b) 


These pole locations correspond to the eigenvalues for the cutoff conditions of surface 
wave modes (energy confined within the dielectric coating). Since the dielectric coating 
which covers the waveguide aperture is located on the ground plane, only the odd TE and 
the even TM surface wave modes can exist. In addition to these poles there is a branch 
point at 0=1; however, it may be noted that 0 = N is not a branch point. 

The contribution of the semicircular contours that occur at each surface wave pole 
may be evaluated by residue theory. The amount of real power confined to the dielectric 
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coating is represented by the surface wave conductance, which may be computed from 
multiplying irj by the sum of the residues of equation (34). This operation gives, for 
the TE surface wave conductance. 



and, for the TM surface wave conductance, 



where the values of for the TE and TM surface wave modes are given by equa- 
tions (35a) and (35b), respectively. 

For the lossy dielectric coating, the values of are complex; thus, no poles 
are located on the real /3-axis and the numerical integration of equation (34) is more 
straightforward. 


PROPERTIES OF THE UNCOATED APERTURE 
Input Admittance 

A partial check on the admittance expression for the dielectric covered aperture 
is that it should revert to the free -space solution when either d = 0 or = e Q . 

For the aperture with no dielectric coating (d * 0 or - Cq), equation (34) 
reduces to 
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which agrees with the expression derived by Mishustin (ref. 1), in which a Green’s func- 
tion approach was used. 

From the phase condition of equations (7), equation (38) may be separated into its 
real and imaginary parts, which results in 



These integrals must be numerically evaluated. 


Asymptotic Expansion of Radiation Fields 

For the aperture with no dielectric coating (d = 0 or q = e Q ), equations (12) 
and (13) reduce to 


Tl(l) = A(1 + r)(x u ) 2 J 1 (x 11 ) 





and 


T 2 (£) = -A(l + DwpoeoJ^Xjj)^^ 

I ■ ^zo 


(40) 


( 41 ) 
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Therefore, for d = 0 or = e 0 , the vector potentials (eqs. (6)) outside the structure 
become 



= aA(l + r)Jj(xu)cos <f> 



(42a) 



-u)M 0 e 0 A(l + HJ^xj^sin 4 > 


^ Jji<€a)£iU P)e^£ 

Jo l ikzo - 


(42b) 


In order to compute the radiation fields from equations (42), the ranges of integra- 
tion must be extended over the entire real axis in mode space. This computation can foe 
made by first noting that as £p — <*>, the Bessel function in the integrand of equations (42) 
can be asymptotically expanded to give 


Jj(|p) 


-3 T 
e * 

iifip 

l/Sir 





(43) 


The substitution of equation (43) into equations (42) results in a sum of two integrals 
with limits over the positive real axis. However, if ij is replaced by in the inte- 
gral containing eS^P, the sums may be combined and the limits become and + 00 , 

A*ni jj. 

Hence, — — and may be written as follows: 




A?" 


-jll 
J 4 


|-=A(l + r)J 1 (x n )^ r cos <# > 




j 1 ’(|a)e 


jr{k zo C Q s sing)] 



/ f 

•a\ 2 l 




-) 

\ji r siif? 


V 

n/ 

J 


(44a) 
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-j— 

A^.-AU + r) 4 



Jl(|a)e 


3 r (k zo c °s 0+^sin 0) 


^k zo ^r sin f? 


(44b) 


where a transformation to spherical coordinates such that z = r cos 8 and p = r sin & 
was made. 

The extension of the integration is required in order to apply stationary phase eval- 
uation of the integral. From Knop (ref, 5), the asymptotic solution of an integral of the 
form 


= g(£)e J 

„ oO 


jrf(£). 


d£ 


(45) 


is given by 


I « g(k Q sin 9)e 


jkpr ^ /-2ir}k 0 


cos 9 


(46) 


When equations (44) are expressed in the form given by equation (46), the far -field 
expansions of the vector potentials become 


A * 111 jk Q r 


r 



cot 9 cos 0Jj f (k o a sin 6) 


/ ,2 



1 - (~~ sin e\ 


l 

u Wl / J 

J 


jk 0 r r 

A m A/1 , ^ VT ^M 0 ^e /sin </> J^a sin Q)\ 
A Z = A(1 + r)j^x n ] — - 


k 0 r 


sin 2 0 


J 


(47) 


It is of interest to note that since the integrands of equations (42) are odd functions 

of i|, f d^Jjdp) . . . may be converted to i \ . d|Hi^(|p) .... Then equa- 
JO 6 '-‘ace 3 ‘ 

tion (16) from Samaddar (ref. 6) with n = 1 may be used to obtain equations (47) directly. 
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For large values of r and sin 6 4 0 ? it is not difficult to show that the radiation 
fields become 


or 


E r m * 0 




TTT A*HI 

E^ 1 * jko —r— sln 8 


„ in _ . . m . 

E 0 * -jwAz sin 0 




r 




in 


J k o r 


“ -JA(1 + ryjifxn) Vi 


k Q a cos 0 cos 0Ji'(k o a sin 0)1 


/k-.a 

1 . _2_ S i n o 
\ x ll 


m J k o r Jsin 4>Jj(k 0 a sin 0)1 

E “ - -jA(l + r) J l( x n )V- *1°, 1 


(48) 


(49) 


The magnitude of the radiation pattern is defined as the square root of the sum of 
the squares of the terms in braces. Although the method and the geometry for the pres- 
ent problem are somewhat different, the above expressions for the radiation fields resem 
ble those obtained previously (ref. 7). The Kirchhoff integral was used several years 
ago (ref. 7; see also, ref. 2) to derive the radiation fields of an open-end waveguide with 
no ground plane. These results can easily be modified to describe the present problem. 


Results 

The input admittance expression (38) was programed on a data processing machine 
in which the Gaussian quadrature method was used to evaluate the integrals. The ratio of 
the diameter to the wavelength 2a/X Q was chosen as the independent variable, and cal- 
culations were performed over the range 0.5859 < 2a/X 0 < 0.9720 which brackets cutoff 
of the TEji mode on one end of the band and that of the TE 21 mode on the other end. 
Although cutoff of the TMqj mode occurs at 2a/X 0 * 0.7654, it is believed that the aper- 
ture discontinuity should not induce a reflected TMqj component, provided a TMqi mode 
is not incident upon the aperture. 
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Results of the computations are shown in figure 2. The input admittance is initially 
large due to the small value of the normalizing guide admittance near cutoff of the TE^ 
mode. However, as 2a A 0 increases toward cutoff of the TE£ j mode, the normalized 
conductance approaches unity, and the inductive susceptance approaches zero. In other 
words, the antenna becomes more nearly matched to free space as the frequency 
increases. 



Figure 2r Input admittance of circular waveguide radiating into free space. 
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The aperture admittance of the circular waveguide, defined as 



is plotted in figure 3. The calculated values show excellent agreement with Mishustin 
(ref. 1). Although Mishustin plots his susceptance on the positive axis, the present cal- 
culations indicate that the susceptance actually goes inductive at the lower end of the 



Figure 3.- Admittance of circular and square waveguide apertures of equal areas 
radiating into free space. 
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Figure Ar Radiation patterns in 0 = tt/ 2 plane for circular and square 
waveguide apertures of equal areas* 
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Figure 5.- Radiation patterns in 0-0 plane for circular and square 
waveguide apertures of equal areas. 




band as verified by preliminary experimental results. The corresponding values of 
yap/yo 316 given lor a square aperture of equal area, excited in the TEqi rectangular 
waveguide mode (ref. 8). The latter is plotted as a function of the guide dimension in 
wavelengths b/X Q where b = ai/?. The conductance of both waveguides is virtually 
identical, and the suseeptanee of both is small over the range considered. This simi- 
larity in the admittance is somewhat surprising in view of the difference in geometry 
and field distribution. 

The normalized radiation patterns in the two principal planes are given in figures 4 
and 5. They, too, are similar to those of the square waveguide (ref. 9), which are indi- 
cated by the broken lines. 


CONCLUSIONS 

The input admittance and radiation patterns of a circular waveguide, opening onto 
a flat ground plane, have been derived by use of the integral transform methods. The 
aperture admittance and patterns for no dielectric coating were computed and are quite 
similar to those of a square waveguide of equal area, excited in the TEqj rectangular 
waveguide mode. 

It is of practical importance to note that, for possible use as an antenna, the open- 
end circular waveguide is well matched to free space above midband of the waveguide. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., April 15, 1968, 

125 - 22 - 02 - 02 - 23 . 


APPENDIX 


FIELD EQUATIONS 


Expansion of the fields in equations (1), (2), and (3) into cylindrical coordinates 
gives, for z £ 0, 


For 0 g z £ d, 


/.♦IN 

p I_ I 8 [ Az 
P ~ P 8 tfAe 0 


A 


w I 3 /"Z 

H -5 ptv 


E z a = 0 


Hn - - r 


i 2 /a 


♦IN 




j^Mo ®P ®z\e 0 

1 1 a 2 fa 1 


3 W P 0 P € o 

I _ 1 1 8 L 8 PlH . 1 1 a 2 

'Z 


-LI0 _82JA£ 

jWju 0 p 8p£ 8 P ye 0 /[ jwp 0 ^2 i*2\e 0 


H 




E 


/ A *H 

n _ 1 a [A g 


‘p ' P a 0 1 ej J j W M 0 ei a p a z 


;K n ) 


A 


e/ . - 


l a 1 


Bp \ € i j i w Po e i p ®0®2 


■K n ) 


r n _ 
^Z - 


H 


n 


j w M 0 e l p ap l 
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Substitution of the vector potentials from equations (4), (5), and (6) yields, for the 
vector components of the fields for z % 0, 
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